Multiple Quantum Phase Transitions at the superconducting LaTiOa/SrTiOs interface 
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We study the magnetic field driven Quantum Phase Transition (QPT) in electrostatically gated 
superconducting LaTiOa/SrTiOa interfaces. Through finite size scaling analysis, we show that it 
belongs to the (2+l)D XY model universality class. The system can be described as a disordered 
array of superconducting islands coupled by a two dimensional electron gas (2DEG). Depending on 
the 2DEG conductance tuned by the gate voltage, the QPT is single (corresponding to the long range 
phase coherence in the whole array) or double (one related to local phase coherence, the other one 
to the array) . By retrieving the coherence length critical exponent u, we show that the QPT can be 
"clean" or "dirty" according to the Harris criteria, depending on whether the phase coherence length 
is smaller or larger than the island size. The overall behaviour is well described by a theoretical 
approach of Spivak et al, in the framework of the fermionic scenario of 2D superconducting QPT. 

Two-dimensional electron gas at the interface between insulating oxides (0-2DEG) have raised considerable 
interest [TH3]- Indeed, they display very high mobilities suitable for applications [5], but also a rich phase diagram, with 
different quantum ground states such as magnetism [SHE] or superconductivity[5J [TU] for example. Moreover, these 
properties can be finely tuned using a gate voltage|TTJ[T2]. Therefore, 0-2DEG appear as interesting systems to study 
Quantum Phases Transitions (QPT) that occur between different quantum states, when a parameter in the Hamil- 
tonian crosses over a critical value [13]. The critical behaviour of the observables belong to universality classes which 
depend on general properties of the system such as its dimensionality or its symmetries, and not on the microscopic 
details. The associated critical exponents obey specific rules, among which the so-called "Harris criteria" . The latter 
stipulates that the correlation length exponent v must satisfy v >2/d for "dirty" disordered systems and v <2/d for 
"clean" ones, where d is the spatial dimensionality [TJ. 

An important example of QPT is the transition from a superconducting to an insulating state in two dimensions, 
which has been a matter of debate for a long time|15j. Numerous experiments with contrasted results 15, have 
been performed, and a large variety of critical exponents have been found. The nature of the non-superconducting 
state and the role of the disorder are still unclear. The possibility of observing multiple QPT has been raised 
recently [TBl - fI5] . Here we show that a perpendicular magnetic field applied on a superconducting 0-2DEG drives the 
system towards a weakly localized metal. We evidence two QPT corresponding to the "dirty" and "clean" limits of 
the Harris criteria, which can be controlled by the gate voltage. More precisely, we have shown that the QPT in 
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superconducting 0-2DEG belongs to the (2+l)D XY universality class|T5J, and follows a fermionic scenario [TOl 120) . 
Finite size scaling analysis (FSS) reveals that depending on the gate voltage and the temperature, the product 
of the critical exponents zv [z is the dynamical exponent) is zv ~ 2/3 in the clean regime in agreement with 
calculations [2 1, and greater than one otherwise, as expected from the Harris criteria. We argue that this multiple 
QPT can be understood in the frame of a recent model developed by Spivak et al.^\ to describe the QPT in a 
disordered array of superconducting puddles coupled by a 2DEG.This result demonstrates that the superconducting 
0-2DEG tunability is a powerful tool to study fundamental properties of superconductors in low dimensions. We 
anticipate that further insight within QPT involving two dimensional superconductors will emerge. For example, the 
role of disorder and inhomogeneitv [T71 122) . Coulomb repulsion, charge fluctuations and screening on the QPT [251 - E5] 
can be explored in more details. 

The LaTiOs/SrTiOa epitaxial interface exhibits a high mobility 0-2DEG which displays superconductivity |10| . 
We recently showed that Tc can be tuned by electrostatic gating from its maximum value of ^200 mK to zero, and 
that superconductivity coincides with the presence of Highly Mobile Carriers (HMC) at the edge of the quantum 
well formed by the 0-2DEG at the interface[T2. It is therefore possible to prepare the system with a given T^. 
by controlling the gate voltage Vq (inset Fig. [T^), and to study how the superconducting state is destroyed by a 
perpendicular magnetic field in this situation. Samples are grown by pulsed laser deposition of 15 unit cells of LaTiOa 
on Ti02-terminated (001)SrTiO3 substrate (see Biscaras et al. 1Q\ for details). A metallic back gate is evaporated at 
the rear of the 500 /im thick SrTiOa substrate and connected to a voltage source (Vq). Standard four-probe resistance 
measurements are made with a current sufficiently low to avoid any heating of the electron at the lowest temperature. 
The polarization scheme described previously [26J is applied to insure a reversible behaviour of the superconducting 
0-2DEG. The normal state displays the logarithmic temperature dependence of the conductivity, characteristic of 
weak localization in two dimensions [TP] . Hence, the magnetic field turns the 2D superconductor into a 2D weakly 
localizing metal as shown in Fig. [T^ for Vq =+80 V, where the resistance per square Rs is plotted as a function of 
temperature T for different perpendicular magnetic fields B. A close-up view of the data reveals a critical field 
which separates the two regimes, and for which Rs is constant (Fig. [T]d). The same set of data plotted as Rs versus 
B for different temperatures in Fig. [2^ displays a crossing point (i?x =372.4 il/D, Bx =0.185 T) where the resistivity 
does not depend on temperature. This is a first signature of a continuous QPT. 

In such zero temperature transitions, the ground state of the Hamiltonian is changed by an external parameter 
(for now magnetic field). Glose to the transition the correlation length ^ in the space dimensions and the dynamical 
correlation length ^t- in the imaginary time dimension of the quantum fluctuations diverge with a power law dependence 
of the distance from the transition 5 — \B — B^ | OIH]. At T=0 the correlation length exponent v defined as ^ cx l^j^"^ 
and the dynamical scaling exponent z defined as oc f ^ are believed to be independent of the microscopic details 
of the transition and only depend on few properties of the system, such as the dimensionality and the range of 
the interactions, which define universality classes for the QPT[T3j. The effective dimensionality of the system is 
d + z where d is the spatial dimensionality. At finite temperature, the imaginary time dimension is limited by the 
temperature fiuctuations so that the dimensionality of the system is d + z at T—Q only, and d at finite temperature|15) . 
More precisely, the finite temperature limits the size of the temporal direction by the thermal cut-off = h/k^T 
which is now an upper bound for the dynamical correlation length near the critical point. It follows that in the 
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spatial dimensions the quantum fluctuations lose phase coherence over a temperature dependent dephasing length 
cx 13]. This leads to a so-called finite size scaling (FSS) of the observables of the system. For instance, 

the resistance takes the form : 



where is a universal function with F(0)=1 ^57]. The critical exponents can then be retrieved by a scaling 
procedure [T5] as follows. The resistance is rewritten as R{6,t) — RxF{\S\t) with t an unknown parameter that de- 
pends only on T. The parameter t is then found at each temperature T by optimizing the collapse around the critical 
point between the curve R{6,t(T)) at temperature T and the curve R{S,t{To)) at the lowest temperature considered 
To, with i(To)=l. The dependence of t with temperature should be a power law of the form t = {T/To)"'^/^'' in order 
to have a physical sense, thus giving the critical exponent product zi>. The interest of this procedure is to perform 
the scaling without knowing the critical exponent beforehand. The result of this procedure applied to the data dis- 
played in Fig.[I^ shows that data collapse onto a single (bi- valued) curve in Fig.[2]D, and yields zv—O.dG (inset Fig.[2]D). 

In the literature, zi>—2/5 has been observed for perpendicular magnetic field driven transitions in conventional 
2D disordered superconductors such as a-NbSi[28l [23 or a-Bismuth|15 1 [30 ] . This value corresponds to the "clean" 
(2+l)D XY universality class according to numerical simulations |2 11 [5T1 132 j . where the extra dimension refers to the 
imaginary time in the quantum transition[13 . In that case, long range superconducting correlations are destroyed by 
quantum phase fluctuations. In general, the dynamical exponent z is found to be one, corresponding to long range 
Coulomb interaction between chargesfT^ I27j. as it has been measured in a-MoGe for instance [20]. In the absence 
of specific screening or dissipation mechanism 33J , there is no reason to invoke short range interactions which would 
set z 1(31]. Therefore, if z = 1, then i/ = 2/3, which means that ly < 2/d, with the spatial dimension d = 2. 
According to the very general Harris criteria[T4], the system is in the "clean" limit at the relevant scale for the 
transition, namely the dephasing length L,^. In 2D disordered systems, i.e. in the "dirty" regime, v is expected to 
increase beyond 1(31] [32], following the Harris criteria, as observed in a-MoGe|^, InOx[2S] or ultra-thin High Tc 
superconductors |36j for instance. 

We now focus on the same Rs vs T curves below 0.1 K (Vg =4-80 V). A blow-up in Fig. [T]: evidences another 
critical field Bq for which the resistance is constant in a restricted range of temperature and ^ changes sign. 
Reported on Fig. Rs vs B curves display a crossing point (5^=0.235 T, i?c=376.6 ri/D), and the FSS analysis 



shows a good collapse of the curves, leading to a critical exponent product zv^l.b (Fig. [2p). This is clearly greater 

than one. The system is in the dirty limit of the Harris criteria at the scale L$. As a conclusion for Vg =-I-80 V, the 
system undergoes a "clean" QPT {zv—2/2i) whose critical behaviour can be seen between 120 mK and 220 mK, and 
a "dirty" QPT {zv—'i/2) at lower temperature. 

When the gate voltage is tuned to lower or even negative values, the system is driven toward a more resistive state, 
and Tc decreases (see inset Fig.[T^). In Fig. [3^, Rs is displayed as a function of the temperature for different magnetic 
fields and for a gate voltage Vq =-15 V. A plateau is clearly seen down to the lowest temperature for a critical field 
i3x=0.1 T (i?x=2176 r^/n), as confirmed by the presence of a crossing point (inset of Fig. [3^). The FSS analysis 
reveals that two regimes take place as shown in Fig. ^p. Indeed, two distinct zv values can be extracted : at "high" 
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temperature (70 mK to 120 mK) z;/=2/3 and at the lowest temperatures zv=?>l2. The corresponding collapses of the 
data depicted in Fig. |3]:,d confirm the existence of the QPT. 

We made the same analysis for all the gate voltages between Vq =-40 V and Vq =+100 V. The results are 
presented in Fig. |4j For Vq >+10 V, two distinct critical fields Sx and Bq [Bq > Bx ) can be found (region II 
in Fig. [4^), whereas for Vq <+10 V, they merge into a unique value (region I). Two important observations need 
to be made : (i) Bx increases with Vq and then saturates when Bx 7^ Be to a value B^- (ii) -Be matches Tc 
defined as the temperature where the normal state resistance drops of 10% [12]. Bq is therefore the critical field 
which fully destroys superconductivity in the system. The critical exponent product zv extracted from the FSS 
around both critical fields are reported as a function of Vq in Fig. |4]d. zi' corresponding to Bx is constant and 
equals to 2/3 : the QPT is in the clean limit. For Bq, zv changes with Vq, with rather important error bars 
due to the finite scaling range in temperature, but is always larger than 1, indicating that the QPT is in the dirty limit. 

To account for these observations, we propose the following scenario based on the XY model where superconductivity 
is destroyed by phase fluctuations in a two-dimensional superconductor. We suppose that the system consists in 
superconducting islands coupled by non-superconducting metallic regions (see Fig. [S] for a sketch). Indeed, as shown 
previously 12: , the superconducting 0-2DEG is made of two types of carriers, a few Highly Mobile Carriers (HMC) 
and a majority of Low Mobile Carriers (LMC). We argued that the presence of HMC triggers superconductivity in 
the system. The density of HMC[T2] (in the 10^^ cm^^ range) and its evolution with the gate voltage is very similar 
to the supcrfluid density directly measured by Bert et [37] in LaAlOa/SrTiOs interfaces. We therefore assume that 
HMC form the supcrfluid with intrinsic inhomogcncity due to its very low average density and the associated density 
fluctuations, and that LMC form the metal which provides long range coupling. The system is therefore described 
as a disordered array of superconducting puddles of size Ld coupled by a metallic 2DEG. Such a situation has been 
studied by Spivak, Oreto and Kivelson (S0K)[T7]. Superconducting phase coherence is governed by the transport 
properties of the metallic part through the proximity effect, and depends on the 2DEG conductance G2deg- If the 
coupling is strong enough (high G2deg)i puddles can develop full local superconductivity and two phase transitions 
can be observed corresponding respectively to the puddles themselves and to the disordered array of puddles. In the 
opposite case (low G2deg)j decoupled puddles are always in a fluctuating regime, and only the full array transition 
can be seen at low temperature. These two regimes can be observed in LaTiOa/SrTiOa interfaces, and the transition 
from one to the other can be controlled by the gate voltage. Indeed, the conductance Gs — l/Rs increases with 
Vq (inset Fig. [T^), and so does the coupling. Region 1 is therefore the low coupling regime with a single transition 
{Bx—Bq), whereas region II refers to the high coupling one with two transitions {Bq > Bx)- In the later case, Bx is 
the critical fleld for the puddle transition and Bq the one for the whole array : as a consequence Tc scales with Bq 
as observed experimentally. 

The SOK model develops a full analysis of the long range coupling between puddles. They introduce the concept 
of "optimal puddle" [17] to account for the statistical distribution of puddles sizes, and show that the critical 
magnetic field Bq scales with the coupling parameter G2deg- At low conductance, Bq cx G2DEG, whereas as high 
conductance, Bq cx G^I^^q. To test the model in more details, we plotted the critical field Bq as a function of 
the conductance G5 on a log-scale in Fig. |4j;. The data are in good agreement with the theory. Not only the two 
regimes can be clearly identified, but also the values of the slopes correspond to the calculated one. This is a strong 
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indication that the SOK model is a good representation of the physics involved in these experiments. 

On this basis, we can now analyse the QPT, remembering that the thermal dephasing length L$ ^ y^i/^ (-^ = 1 
here). A cartoon illustrating the situation is shown in Fig. [sj Let us first focus on region II as defined in Fig. |4] 
At high temperature, the phase correlation length is smaller than the superconducting puddles size Lfi, and the 
intra-puddle physics dominates. The critical field Bx corresponds to the dephasing field Bd on a puddle of size Ld 
{Bd ^ *&o/id where $o is the flux quantum) [17l [18], and does not depend on the microscopic parameters of the 
system tuned by Vq. This is why it is constant as a function of Vq (B^ = Bd) as shown in Fig. |4^. From the 
value of Bd, we can estimate Ld to be of the order of 100 nm. Since L$ < Ld, the system is in the clean limit, and 
zv — 2/3 < 1 (Fig. [Ija) as expected from the Harris criteria. When lowering the temperature, L^ crosses Ld at Td, 
and the whole array undergoes the transition with a critical field Bq- In that case, phase fluctuations extend over 
wide disordered regions and the dirty (2+l)D XY model applies (zv > 1 as shown in Fig. |4|d). In region I, only the 
array displays a transition with a single critical field Bq. However, as the temperature is lowered, L$ also crosses 
Ld, and we therefore observe a transition from a clean {zi^ = 2/3 < 1) to a dirty {zv > 1) (2+l)D XY limit. At 
this point, it is worthwhile mentioning that, strictly speaking, the transition at Bq is the only true QPT, since the 
scaling holds down to the lowest temperature. However, at sufficiently high temperature, the transition at Bx displays 
fluctuations corresponding also to a real QPT which is never reached. Indeed, the diverging coherence length crosses 
Ld, which acts as a cut-off length for the intra-puddle physics. Bx is therefore a crossover fLeld[T7], but we consider 
both transitions as QPT in view of the quantum nature of the fluctuations. 

Such an analysis reveals that the tunable LaTiOa/SrTiOs epitaxial interface appears as a unique system to 
study superconducting QPT in 2D, which is well described by the (2-|-l)D XY model as expected. The role of 
disorder on the critical exponents is clearly evidenced, together with the possibility of observing multiple QPT as 
previously proposed theoretically [TBHT5] . These models have been developed in the framework of the so-called 
"fcrmionic scenario" [13 [2D] , where superconductivity disappears because Cooper pairs are destroyed, as opposed to 
the "bosonic scenario" '27], where they localize and form an insulator. In the later case, the QPT takes place for 
critical resistances Rc close to the quantum resistance Rq — h/Ae^ ~ 6.5 k57/n, whereas in the former one it can 
occurs at much lower resistances 1 1 6IH8] . Steiner et a/. [35J compiled experimental results in the literature and clearly 
evidenced the two behaviours : the QPT separates a superconducting phase from a weakly localizing metal in the 
less disordered systems [Rq < Rq), and from an insulator in strongly disordered materials (i?c ^ Rq)- Our data 
(0.04 < (Rq/Rq) < 0.4) fully agree with this picture, and the non-superconducting phase is indeed always a weakly 
localizing metal. 

As a conclusion, we showed that the superconducting 0-2DEG at the LaTiOa/SrTiOa interface undergoes a QPT 
from a superconductor to a weakly localizing metal (fcrmionic scenario) upon applying a perpendicular magnetic 
field, driven by phase fluctuations and well described by the (2+l)D XY model, as expected. By tuning the gate 
voltage, it is possible to explore the "clean" and "dirty" regimes according to the Harris criteria, with a critical 
exponents product zv — 2/3 in the former case, in agreement with previous evaluations |11| [ 55] . and greater than 
one otherwise. The system is well described by a disordered array of superconducting puddles coupled by a 2DEG, 
which can exhibit two QPT, one related to "regional" or local ordering, and another one corresponding to long range 
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phase coherence, as proposed theoreticallv |16f[T5] . The key parameter, that is the couphng constant (the 2DEG 
conductance), can be tuned at will to explore the phase diagram of the system. This is important in the context 
of recent studies of strongly disordered 2D superconductors where intrinsic inhomogeneities appear at mesoscopic 
scales pn 1391 140] . with coexistence of superconducting and non-superconducting regions. Recent work on artificial 
ordered metallic networks addresses this issue|41j. but does not reach the insulating state. Our study opens the 
way for exploring the physics of disordered superconductors, and beyond, the more general problem of phase co- 
herence in multiscale systems such as strongly correlated materials which are phase separated |42j or spin-textured [43] . 
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Figure 1: Superconductor-insulator transition induced by a magnetic field, a, Sheet resistance Rs as a function 
of temperature for different magnetic fields from to 0.3 T. Inset : Tc and Gs~l/Rs as a function of the gate voltage 
Vg- b, Zoom on the same data showing the characteristic magnetic field Bx, for which Rs is constant between 0.12 K and 
0.22 K. c, Zoom on the same data showing the characteristic magnetic field Be, which separates the two regimes at the lowest 
temperatures. 
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Figure 2: Finite size scaling analysis for Vg =+80 V. a, Sheet resistance Rs as a function of magnetic field B for different 
temperatures from 0.1 to 0.2 K. The crossing point is (Bx =0.185 T, 7?x =372.4 f2/n).b, Finite size scaling plot Rs/Rx as 
a function of \B — Bx\t (see text for the definition of t). Inset : temperature behaviour of the scaling parameter t. The 
power law fit gives zv=Q.&Q. c, Sheet resistance Rs as a function of magnetic field B for different temperatures from 0.04 to 
0.07 K (Vg =+80 V). The crossing point is (Bc=0.235 T, i?c=376.6 fi/D). d, Finite size scaling plot Rs/Rc as a function of 
\B — Bc\t. Inset : temperature behaviour of the scaling parameter t. The power law fit gives zi'=1.5 
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Figure 3: Finite size scaling for Vg =-15 V. a, Sheet resistance Rs as a function of temperature for different magnetic 
fields from to 0.2 T. Inset : corresponding Rs as a function of the magnetic field for different temperatures from 0.07 K to 
0.12 K. The crossing point is (i3x=0.1 T, _Rx=2176 f^/D) . b, Temperature behaviour of the scaling parameter t : two distinct 
slopes are evidenced, 0.64 at high temperature and 1.4 at low temperature, c, Finite size scaling plot Rs/Rx as a function of 
\B — Bx \t corresponding to the "high temperature" regime, with zh'=0.64 . d, Finite size scaling plot Rs/ Rc as a function of 
\B — Bc\t corresponding to the "low temperature" regime, with zu—lA 
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Figure 4: The two quantum phase transitions, a, Bx , Be (left scale) and Tc (right scale) as a function of Vg- The dashed 
line corresponds to Bd (see text). Regions I and II refers to the low, and respectively high, coupling regimes (see text), b, zu 
as a function of Vg for the two transitions (when Bx=Bc, Be has been assigned to the "low temperature" regime, and Bx to 
the "high temperature" one), zu is mainly constant with zv=Q.6Q for the "high temperature" transition. It is between 1.5 and 
2.3 for the dirty one. The dashed lines correspond to zv=2/?i and zi/=3/2. c, Bq as a function of Gs on a log scale. Dashed 
lines have slopes 1 and 1/4 respectively as expected from the SOK model. 



13 




Figure 5: Sketch of the role of mesoscopic disorder on the QPT. The dephasing length diverges with temperature, 
and reaches the size of the superconducting puddles Ld at Td- The insets show a piece of material in the two regimes. At high 
temperature (bottom), Z/$ < Ld and the system is in the clean limit, whereas at low temperature (top), L$ > Ld, and the 
system is in the dirty limit. In this drawing, superconducting puddles (in blue) are coupled through a 2DEG (in yellow). The 
arrows symbolize the local phase of the superconductor 



